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Abstract 



The exponent A that describes the decay of the autocorrelation function A(t) in 

a phase ordering system, A(t) ~ L^^'^~^\ where d is the dimension and L the 

Tjj- ■ characteristic length scale at time t, is calculated exactly for the time-dependent 

^ ■ Ginzburg-Landau equation m. d = 1. We find A = 0.399 383 5. . .. We also show 

c^ \ explicitly that a small bias of positive domains over negative gives a magnetization 

which grows in time as M{t) ~ L'^ and prove that for the Id Ginzburg-Landau 

'^ ! equation, // = A, exemplifying a general result. 
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The field of pliase ordering kinetics lias seen a number of new developments 
in recent years |]I|. In particular the values of the growth exponents z, which de- 
scribe the time-dependence of the characteristic scale L{t) via L ~ t^'^, are known 
exactly for most models with purely dissipative dynamics |]I], 0]. For systems with 
short-range interactions and dynamics which are either nonconserved or obey a local 
conservation law, the exponent z is usually a dimension- independent integer 0. Re- 
cently, however, it has been realized that for nonconserved dynamics the description 
of two-time correlations requires a new exponent, whose dependence on the spatial 
dimension d and on the symmetry of the order parameter is nontrivial [§, ^ . 

The exponent A can be defined in terms of the general two-point correlation func- 
tion C(r; ti, ^2) = (0(x, ti) 0(x-t- r, ^2)), where is the order parameter field. In the 
scaling regime, this is expected to have the scaling form C{r] ti, ^2) = f{r/Li, r/L2), 
where Li, L2 are the characteristic length scales at times ti and ^2 0) 0- Iii the 
limit of well-separated times, L2 ^ Li, one anticipates 0] the power-law form 
C(r;ti,t2) ~ {Li/ L2Y^^ f{r / L2), defining the exponent A. An especially simple 
case is where we take r = 0, and the initial time ti = 0. Then the general form 
reduces to 

A(t) = C(0;0,t)~[W^(t)f"\ (1) 

where ^0 is some fixed length related to the initial conditions. The 'autocorrelation 
function' A{t) has been measured in simulations of 0{n) models for various spatial 
dimensions (i P, |^, and in experiments on twisted nematic liquid crystals films [§, 
and the exponent A deduced. It generally has a nontrivial value. 

There are a few analytical results for A - the nonconserved 0{n) model forn = cxd 
(A = (i/2 [^), and the d = 1 Glauber model (A = [Q), while for nonconserved 
scalar fields in d = 2 Fisher and Huse [^ have conjectured that A = 3/4 exactly. 
In general, however, A appears to be a nontrivial exponent associated with ordering 
dynamics, although it is known to satisfy the bound (in our notation) A < d/2 for 
nonconserved dynamics p, [1^ . 



In this paper we calculate A exactly for a soluble model corresponding to the 
late-time, zero-temperature coarsening dynamics of the time-dependent Ginzburg- 
Landau (TDGL) equation for a scalar field in d = 1. The equation of motion is 
dt(j) = d'^cj) — dV/dcj), where V{(j)) is a symmetric, double-well potential with minima 
at = ±1 (e.g. V^(0) = (1 — (f^Y)- At late times, when the mean separation L of 
domain walls is large compared to their intrinsic width ^ (= [V^"(l)]~^/^), the walls 
only interact weakly, through the exponential tails of the wall profile function. Then 
the dynamics is very simple [0, |1^, |T^. The closest pair of walls move together 
and annihilate, while the other walls hardly move at all, and the system coarsens by 
successively eliminating the smallest domains. It is found that the distribution of 
domain sizes / approaches a scaling form, P[l) = L^^f{l/L). The scaling function 



/(x) can be exactly calculated [jTT], |T2|, ^ 



In an earlier work [O], we have shown that there is a nontrivial exponent asso- 



dated with the fraction of the hne that has never been traversed by a domain wall 
(i.e. the fraction of the line where the order parameter cj) has never changed its sign 
ig). This fraction decays as L-^'^-^\ with p = 0.824 924 12 . . .. Here we show that 
the approach developed in [|1^] can be generalised to calculate A for this model. The 



result is A = 0.399 383 5 . . .. A recent simulation of the same model |T5| gave the 
estimate A = 0.43 ±0.01, which, we think, is in reasonable agreement with our exact 
result, given that the extrapolation to large L was not straightforward. 

The exponent A can also be obtained from the rate at which a small initial bias in 
the order parameter grows with time |]l6l, (0) ~ L'*'. We demonstrate this explicitly 



within the present model in the second part of this work. 

The calculation of the autocorrelation exponent A follows closely that presented 
in reference [0. One starts with random intervals on the line. Each interval / 



is characterised by its length 1{I) and by its overlap q{I) with its initial condition 
(initially q{I) = 1{I) for all /). At each iteration step, the smallest interval /min 
is removed (i.e. the field is replaced by —0 in this interval). So three intervals 
(the smallest interval /min and its two neighbors Ji and I2) are replaced by a single 
interval /. The length and the overlap of the new interval / are given by 

/(/) = /(/i)+/(/mi„)+/(/2) , (2) 

q{I) = q(I,) + q(I,) - q(I^.J . (3) 

Then the average length L of domains and the autocorrelation function A are given 
by 

A = E^(^)/E^(^)- (4) 

where the sums are over all the intervals / present in the system. 

The argument showing that no correlations develop if none are present initially 
was given earlier [|l^ and the calculation is then very similar to that for the eval- 
uation of the exponent (3. We take, for simplicity, the lengths of the intervals to 
be integers and io to be the minimal length in the system. We also assume that 
the total number A^ of intervals is very large. We call rij the number of intervals of 
length i and qi the average overlap of the intervals of length i. At the beginning, 

qi = i- 

We denote with a prime the values of these quantities after all the nj^ intervals 
of length iQ have been eliminated, so that the minimal length has become zq + 1. 
Then the time evolution is given by (compare equation (2) of |T^) 

N' = N-2n,, 

I \^ _Z ""i--3-i'0 

J=«0 

+ "-io E ]^ ' ^ '" (^i + ^i~J-io -lio) ■ (5) 

j=io 
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This is only valid under the condition that nj^ -C A^ which is indeed valid when ig 
becomes large and as long as the system consists of a large number of intervals. 

We assume that after many iterations, i.e. when ig becomes large, a scaling limit 
is reached where 
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where A is the exponent we want to calculate (Q). Because io is so large, we can 
consider x = i/iq as a continuous variable. This gives 



N' 
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«o 



N 



n. 



^0 + 1 

4 = N'{^o + lr-'g(-!— 

\io + 1 



fix) - -f{l)f{x) - -f{x) - -xf\x) 

«0 ^0 «0 
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Inserting these expressions in the time evolution equations (|^) gives 
df 



(7) 
(8) 






dio 
dg 

dio 



fix) + xf'ix) + 9{x - 3)/(l) / dy f{y)f{x - y - 1) 



x-2 



(1 - X)g{x) + xg'{x) + 2d{x - 3)/(l) / dy g{y)f{x - y - 1) 



x-2 



x-2 



g{l)e{x-3) / dyf{y)f{x-y-l). 



(9) 



In (H), both rii and n^gj are functions of x = i/io and of io, and the partial derivatives 
in d^) mean the derivative with respect to io, keeping x fixed. Demanding that the 
system is self-similar, i.e. that the functions f{x) and g{x) do not change with 
time (i.e. replacing the left-hand sides of (^ by zero), one finds that the Laplace 
transforms 



0(P) 
V'(P) 



-px 



-px 



f{x) dx 
g{x) dx 



(10) 



satisfy the following equations (where primes now indicate derivatives) 



-/(l)e-^-p0'(p) + /(l)e-V(p) = (11) 
-X'il;{p)-g{l)e-^-pij'{p) + 2f{l)e~P(P{p)ilj{p)-g{l)e-y\p) = 0. (12) 



Defining the function h{p) by 



hip) = 2/(1) 



-dt , 



(13) 



the solutions of the above equations are 
(j){p) = tanh[h{p)/2] 



ip{p) 



PC 

Jp 



0'(?)) 



1 - 02 (p) q 



x~i 



'dq . 



(14) 
(15) 



1 — 02(g) p^ 

The constants of integration imphed by these forms were fixed by the requirement 
that both and ip decay fast enough for large p, as is clear from the definitions 



(p!0|). So far the parameters /(I), g{l) and A are arbitrary. We shall see that they 
are fixed by physical considerations. 

Eq. ( [I^ ) for 0, which determines the domain size distribution, is of course iden- 
tical to that obtained in previous work |Tl], |12|, 0. Eq. (|1^) for tjj can be rewritten 
in the more convenient form 



ij{p) = 2g{l) 



Ji{q) 



+ e 



-h{q) 



q 



A-l 



,p (M{p) + 2 + e-'^^p) p^ 
It is helpful to introduce the expansion 



'dq. 
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/ dq = -\ogp--f-J2 

-'P ^ n=l 



[-pT 



n n\ 



(16) 



(17) 



where 7 = — /o°° dt e~* logt = .577215 6... is Euler's constant. 

From the small-p expansion of (|14[), it is easy to show that, provided the first 
moment of the domain size distribution exists, one must have /(I) = 1/2 |Tl], |T2|, |13 



From now on, we will consider only this case (see |]T^ for the discussion of cases where 
the stationary distribution has long tails). Defining the function r{p) by 



r{p) = h{p) + \ogp 

jp 

one obtains, using (|16[), 

V'(P) = 2^(1 



-dq + logp 



or{q) 



+ q^e 



-r(q) 



,\-2 



^dq 



;i8) 



(19) 



p e''(p) + 2p + p2e-Kp) pA-i 

Now r{p) can be expanded in powers of p, using (|T^), and so this last form makes it 
easier to analyse the singular behavior of ip{p) aX p = 0. One finds that, for small p, 

i:{p) = A + Bp'-^[l + 0{p)] , 



(20) 



where A = 25f(l)/(l - A) and 



B 



2g{l)e 



-r(0) 



/ ^, \ {r'{q)-l)e''-''^dq+ / gVV 
JO 1 — A JO 



-r(q) 



dq 



(21) 



2^(l)e^(l-A)-W g^"V 
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q 
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+ q\l - A)e-"WJ dq . 



Now compare ( pO|) with a direct expansion of ([lOD , namely ip{p) = J^ dxg{x){l — 
px + 0(p2)). If the function g{x) is to have a finite first moment then we must have 
i? = in (Eol). This condition determines A as 



A = .399 383 5 



(22) 



From numerical simulations of the same model, Majumdar and Huse [^ found the 
power-law decay A{t) ~ L~^, with A = 0.57 ± 0.01, corresponding to A = d — A = 
0.43 ±.01. There were, however, large corrections to scaling in their numerical data, 
which we think are the origin of the disagreement between their numerical estimate 
and our exact result. 

As in []TB[, one can show that B ^ would correspond to a power law decay in 



g{x) and that such a power law cannot be produced if it is not present in the initial 
condition. Note that g{l) cannot be determined as one can always multiply all the 
qi by a constant without changing our results. 

For the remainder of this paper we will look at a related quantity, the growth of 
an initially small bias in the order parameter, and show that the bias grows as L^ 
as the system coarsens (while the bias remains small). Furthermore, we will show 



explicitly that /x = A for this model, exemplifying a general result |]T6 . 

Consider a sequence of positive and negative domains on a line. We call rii 
{rrii) the number of positive (negative) domains of length i. The total number N 
of positive domains is of course equal to the total number of negative domains, 
N = J2i^i = J2i^T^i- When the domains of size ^o are removed, the new values of rii, 
rrii and A^ are given by 



n' 



-L Ar M" '"»o Z^ 






N' = N-m,- rrii, ■ (23) 

Let us write forms for rii and rrii analogous to the first of equations (|) : 

,, = ^/,(l), ™. = ^/,(l). (24) 



Then one has 



/ N -rij,- rrij, ( i ^ 
«o + 1 V^o + 1' 



which gives, for i^ large (when x = i/io can be treated as a continuous variable), 

, N 



ni = — 



h{x) + - {-/i(l)/i(x) - /2(l)/i(x) - h{x) - xf[{x)} 

«0 



(26) 



and a similar expression for m'i. 

Inserting the forms (|2^ ) into ( |26D gives coupled evolution equations for /i and 

^0^1^ = [h{l)-f2{l)]h{x)+h{x)+xf[{x) 

+e{x - 3)/2(l) f^~' dy /i(i/)/i(x -y-l), (27) 



and a second equation obtained by interchanging the subscripts '1' and '2'. Note 
that the derivatives on the left-hand sides are with respect to the (imphcit) second 
argument io- Introducing the Laplace transforms with respect to the first argument, 

/oo 
Ux)e-P^dx, in = l,2) (28) 

one finds that their evolution is given by 

,0^^ = [/,(!) _ f,il)]i;,ip) - pi,[{p) - /i(l)e-^ + /2(l)e-^^?(p) , (29) 

and a second equation with '1' and '2' interchanged. 

So far this is completely general. The basic idea is to perform a linear stability 
analysis around the 'symmetric' solution '?/'i(p) = ''p2{p) = 0(p)) where 0(p) satisfies 
(p!T|) with /(I) = 1/2, in order to determine the rate at which a small perturbation 
will grow. We therefore take ipi{p) and il'2{p) to have the forms 

^n(p) = 0(p)±ea(p), (30) 

with 

fnil) = l±ea, (31) 

with e small and the + (— ) sign corresponding to n = 1 (n = 2). If the bias repre- 
sented by the terms in e is a relevant perturbation, a{p) will grow under iteration: 
a ~ (^o)'^ with yU > (and similarly, a ~ (^o)^ in (|3TD ). Subtracting from (^9]) its 
counterpart with '1' and '2' interchanged, and putting iod(r{p)/dio = /icr(p), yields 
the eigenvalue equation 

fia = 2a(j) - pa' ~ ae'^ ~ ae^^^^ + e^^^a , (32) 

with solution 

a{p) = a / - -^rrdq ■ 33 

Jp q \pj l-(j)^{q) 

The integration constant was fixed as before by the requirement that a{p) decrease 



as exp(— p)/p for large p, which follows from (^Sf ), (|30| ) and (|3l|). Demanding once 
more that a{p) be regular at p = (so that the first moments of fi{x) and f2{x) 
exist) yields the following equation for fj,: 

j^ rfg[e-V(g) + e-''-20(g)]^^^^ = O. (34) 

Using 0(g) = (e'"'-'^-' — q)/{e^^'^'' + g), which follows from ( [I^) and (p!^), gives the 
condition 

/■oo 

/ rfg [ (e~^ - 1) e''^^) + q^e"' + 1) e-''^^^ Ig^"' = (35) 

JO 

for /i, with solution yU ~ .399 38.... Comparison with ( ^2]) suggests that /i = A. In 
fact, using integration by parts one can show that condition (^) for fi is identical 
to ( piD (with 5 = 0) for A, and so yU = A exactly. 



The result // = A is, in fact, quite general. For TDGL dynamics, it has been 



discussed elsewhere |^. Let us derive it for any kind of dynamics of an Ising model. 
Consider a system of N Ising spins in dimension d. We call P( {S'j(t)} | {S'j(O)} ) the 
probability of findingthe system in the spin configuration {Si{t)} at time t given 
that it was in configuration {^^(O)} at time 0. We assume that the system evolves 
in a zero magnetic field and that the dynamics preserves the it symmetry, namely 

p{ {s,{t) }\ {sm} ) = p{ {-s^m I {-^.(0)} ). 

Suppose that one starts with an initial condition {S'j(O)} chosen completely at 
random, then the correlation {Si{t)Sj{0)) is given by 

{s^{t)s,m = ^ E E sMSM p{ {s^m I {sm} ) • (36) 

{s{t)} {sm 

where I]{5(j)} indicates a sum over the 2^ configurations at time t. 

Suppose on the other hand that one starts with a weakly magnetized initial 
condition, i.e. the initial configuration {S'j(O)} is chosen with probability 

i=l 

when m{Q) is infinitesimal. Then the magnetization m{t) per spin at time t is a 
function of m(0), and to first order in powers of m(0) one has 

Mt) = E T. n{sm]\{sm])Q{{s^m)^^'^^ 

{S(t)} {5(0)} ^^ 

~ m(0) '—j^ . (37) 

Therefore if one assumes that due to some coarsening phenomenon the two-point 
function scales as 

where Rij is the distance between sites i and j, one finds that 

m{t) ^L^m(O) f d'^RfiR) 

which means that the magnetization and the autocorrelation exponents are the same. 

To summarise, we have derived a non-trivial value for the exponent A within an 
exactly soluble model, and shown explicitly that the growth of an initial bias in the 
order parameter is controlled by the same exponent. 

We thank the Isaac Newton Institute, Cambridge, where this work was carried 
out, for its hospitality. 
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